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ABSTRACT 

We report on the optical properties of a sample of 30 BL Lac host galaxies 
in the redshift range 0.03 < z < 0.2, as derived from HST observations. All 
galaxies are fully resolved in the WFPC2 (F702W filter) images, allowing a 
quantitative analysis in two dimensions. Most and possibly all these galaxies 
have characteristics very similar to those of "normal" giant ellipticals. The 
luminosity, ellipticity, isophote twisting and amount of disky or boxy isophotes 
are consistent with those found in non-active ellipticals and in radio galaxies. 
In all cases the BL Lac nucleus is well centered in the main body of its host 
galaxy, a result that argues strongly against the microlensing hypothesis for 
any significant fraction of the population. A search for faint sub-structures 
in the host galaxies has not revealed notable signatures of tidal distortions 
or sub-components (faint disks, bars, X features, etc.), and with only one 
exception, there are no prominent dusty features in the central regions. Instead, 
the BL Lac host galaxies are smooth and unperturbed, suggesting that strong 
external gravitational interactions are not important to ongoing activity. A 
careful examination of the environment around the nucleus, however, shows 
a high incidence of close companion objects, whose nature remains unclear 
pending spectroscopic observations. 

Subject Headings: BL Lacertae objects:general — BL Lacertae objects: host galaxies 
Galaxies: structure — Galaxies: interactions 
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1. Introduction 

The Hubble Space Telescope (HST) snapshot imaging survey of BL Lac objects^ has 
produced a large, homogeneous data set of over one hundred high-resolution images for the 
study of the environments of BL Lac nuclei. The targets were drawn from various samples 
and include sources in a wide range of redshift (z from 0.03 to ~ 1.4) as well as objects 
with unknown z. Results from the survey are given in Scarpa et al. (2000; hereafter, Paper 
I), and Urry et al. (2000; hereafter, Paper II). Previous HST imaging of BL Lacs had been 
done only for a small number of sources (see Falomo et al. 1997; Jannuzi et al. 1997; Urry 
et al. 1999) but demonstrated the capabilities of HST for studying the environs of BL Lacs. 

The HST snapshot images have two main advantages with respect to ground-based 
images (e.g., Falomo 1996; Wurtz, Stocke, & Yee 1996; Falomo & Kotilainen 1999): an 
order of magnitude improved spatial resolution and very good homogeneity of data. On 
the other hand the snapshot images were obtained with relatively short exposure times 
(typically minutes) and therefore are not very deep. This sets limits to the detection of very 
faint extended features, a well-known problem for HST images (see Falomo & Kotilainen 
1999 for a direct comparison of HST and ground-based images of BL Lacs). 

Our first papers on the overall properties of the BL Lac objects derived from HST 
data focused on the determination of the absolute magnitude, scale-length, and type of 
the host galaxies, as well as the luminosity of the nuclear source (see Paper II). This was 
based on accurate decomposition of the azimuthally averaged luminosity radial profiles 
which allowed the study of the sources down to surface brightness jiR ~25 mag arcsec~ 2 . 
Additional morphological properties of the surrounding nebulosity, such as ellipticity, 
isophotal twisting, isophote shape (disky versus boxy), presence of dust, etc., can only be 
derived from a full two-dimensional approach, which is possible for a subset of bright nearby 
host galaxies. 

For elliptical galaxies it has been shown that isophote shapes are related to internal 
kinematics and to the radio and X-ray emission properties (e.g. Bender, Dobereiner, & 
Mollenhoff 1987). Moreover significant distortions of the isophotes and/or twisting may be 
produced by strong galaxy interactions and/or merger events. It is therefore of interest 
to investigate these observable features in the galaxies hosting BL Lac nuclei in order to 
understand the (possible) role of galaxy structure on the nuclear activity. This kind of 
analysis has been performed so far only for a small number of nearby BL Lac objects (e.g., 
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Falomo 1991,1996; Heidt et al. 1999) because the presence of the bright nucleus often 
hinders a fully detailed morphological study. 

In this paper we present the results of a two-dimensional analysis for 30 nearby 
(z ^ 0.2) snapshot BL Lacs, i.e., those with the highest signal-to-noise ratios and the largest 
extent (see Table 1). For 25 of the 110 objects observed in the HST snapshot survey (Paper 
I) the redshift is unknown. Of these 17 are spatially unresolved and are therefore most 
probably at z > 0.2. For the remaining 8 resolved objects with unknown redshift, based 
on the apparent magnitudes of the nebulosities (see Paper I) and assuming the luminosity 
of the host is consistent with that of other snap BL Lacs (see Fig. 1), one expects that at 
most a few can be at z < 0.2. Therefore the objects discussed here represent practically the 
totality of z < 0.2 BL Lacs from our original sample of 110 objects. 

For these low- z objects it is possible with HST to investigate features and structures 
in the host galaxy that are undetectable with ground-based observations. The host galaxy 
can be investigated down to less than 1 kpc from the nucleus (corresponding to 0'.'3 at 
z = 0.2). In particular one can search for slight off-centering of the nucleus with respect to 
the galaxy, determine the presence of sub-components in the host galaxy, investigate the 
detailed shape of the isophotes, and search for dusty features and close companions. 

Because the two broad components in BL Lac spectral energy distributions (SEDs) have 
peak power outputs iyL v ) at wavelengths that increase systematically with luminosity (e.g., 
Ulrich, Maraschi & Urry 1997), BL Lacs found in radio and X-ray surveys have generally 
different SEDs. Radio-selected objects tend to have SEDs peaking at infrared-optical 
wavelengths and in the MeV-GeV gamma-ray band (low-frequency-peaked BL Lacs, or 
LBL) and exhibit luminosities approaching those of quasars. In contrast, X-ray-selected 
BL Lacs have SEDs peaking at UV- X-ray wavelengths and again at TeV energies 
(high-frequency-peaked BL Lacs, or HBL) and are generally less luminous. Of the 30 
sources in the present study, 6 are of LBL type , the rest being HBL (Padovani & Giommi 
1995). In the following analysis we have not addressed whether morphological properties 
are different for the two types of BL Lacs, having found from our earlier analysis that the 
overall host galaxy characteristics do not depend on type (Paper II). 

In § 2 we briefly describe the observations and our data analysis. § 3 gives the results 
on the host galaxies properties and shows the comparison with non-active ellipticals and 
radio galaxies. Finally in § 4 we summarize the results and draw the main conclusions from 
this study. Throughout the paper we used H = 50 km s" 1 kpc" 1 and q = 0. 
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2. Observations and Data Analysis 

The snapshot observations of BL Lacs have been described in detail in Paper I, 
therefore only the main steps are reviewed here. All objects considered here were observed 
with the WFPC2 camera of the HST in the F702W filter. Each final image was obtained 
from the combination of typically three exposures for a total integration time of about 
5-20 minutes (see Table 1 in Paper I for details). The Point Spread Function (PSF) used 
is a combination of a model derived from Tiny Tim package (Krist 1995) plus a diffuse 
exponential tail that accounts for the first order scattered light (see Urry et al. 1999 and 
Paper I for details). Photometric calibration was taken from Holtzman et al. (1995). 

We performed a two-dimensional surface photometry analysis using an interactive 
numerical mapping package (AIAP; Fasano 1994), which produces an isophotal map of the 
image. This map was masked in order to avoid regions contaminated by companion objects, 
diffraction spikes, and any other extraneous features visible in the images. The isophotes 
were then fitted with ellipses down to fiR = 22-23 mag arcsec -2 (see Fig. 2 for an example). 
There are five geometric free parameters per isophote: semi-major axis (a), center position, 
ellipticity (e), and position angle (PA), which allow us to characterize the morphological 
and the photometric properties of the galaxy. 

From this analysis we then derived the average position angles and ellipticity as a 
function of the generalized radius, r = a x (1 — e) 1 / 2 , and investigated their dependence 
on r (isophote twisting and ellipticity profiles). In addition we examined the possible 
displacement of the centers of isophotes and the centering of the nucleus with respect to the 
galaxy. 

3. Results 

3.1. Luminosity Profiles and Two Dimensional Modeling 

For each object we derived a luminosity profile from the two-dimensional isophotal 
analysis and compared with that derived from the azimuthal average reported in Paper I. 
This was possible down to fiR ~ 22-23 mag arcsec -2 since the fainter outer regions are too 
noisy to be investigated with the isophotal map. For all sources we find that the luminosity 
profile derived from the fit with ellipses of the isophotes is always in excellent agreement 
with that extracted from the azimuthal average (see example in Fig. 3). 

We also fitted the isophotes with a two-dimensional galaxy model (either a 
de Vaucouleurs law or an exponential disk) convolved with the PSF plus a point source 
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modeled by a scaled PSF. In all cases the de Vaucouleurs law gives a better fit to the 
data than the exponential disk model, in agreement with Papers I and II. This result is 
supported by the lack of spiral structures seen in the WFPC2 images and by the diskiness of 
the isophotes (see discussion in the following paragraphs). Absolute magnitudes and scale 
lengths of the host galaxies derived from the two-dimensional fitting are also in excellent 
agreement with the results of one-dimensional profile fitting (Paper I), and are consistent 
with the results from our deep HST images for a small number of BL Lac objects (Falomo 



Since the ellipticity e is in general dependent on the distance from the galaxy center, 
we report its value at the effective radius. However, because of the relatively high surface 
brightness limit, in some objects the ellipticity at the effective radius is not measurable and 
in these cases we took the maximum ellipticity. 

The values of e are reported in Table 1. The error on observed e for isophotes far from 
the nuclear regions is typically 0.02 (see Fasano & Bonoli 1990 for discussion). We show in 
Figure 4 the distribution of ellipticity for our 30 low-redshift BL Lacs compared with that 
of 200 nearby normal (radio-quiet) ellipticals (Fasano & Vio 1991) and of a sample of 79 
low-redshift radio galaxies (Govoni et al. 2000). Both normal ellipticals and radio galaxies 
exhibit the same ellipticity distribution: < e >eu— 0.22 ± 0.13 and < e >rg= 0.21 ± 0.12. 
The average ellipticity for BL Lac hosts, < e >bll — 0.16 ± 0.09, is slightly smaller than 
that of non-radio ellipticals and radio galaxies but is still consistent with being drawn from 
the same population (Prs ~ 0.1). 

3.3. Isophotal Displacement and Nuclear De-Centering 

Displacement of isophotes with respect to the center may represent either a global 
distortion of the host galaxy, due to recent tidal interactions with close companions (e.g., 
Aguilar & White 1986), or a consequence of gravitational microlensing (Ostriker and Vietri 
1985). In the first case one would expect non-concentric isophotes to exhibit an asymmetry 
of distribution. To quantify this asymmetry we have computed the dimensionless parameter 



5 = y (X c — X Q ) 2 + (Y c — Y ) 2 /r, where X c , Y c are the centers of isophotal ellipses, X Q1 Y Q is 
the location of the nucleus, and r is the radial distance to the particular elliptical isophote 
(defined earlier). Since 5 may slightly change with r we took the value at the effective 



et al. 1997). 



3.2. Ellipticity 




- 6- 



radius R e of the galaxy or, if R e is beyond our measurements, we took the average value 
excluding the values at r < 1 arcsec. 

In Figure 5 we show the distribution of 5 derived from our sample of low-redshift BL 
Lacs compared with that of radio galaxies from Govoni et al. (2000). The average value 
for our sample is: < 5 >bll= 0.027 ± 0.01 (median 5 is 0.030), marginally consistent 
with zero but suggestive of small distortions in the galaxy isophotes. This is similar to 
the value reported for low z radio galaxies, < 5 >rg= 0.03 ± 0.04 (Govoni et al. 2000), 
using a homogeneous procedure, while for 40 normal ellipticals observed by Sparks et al. 
(1991) from the ground, Colina & de Juan (1995) give < 5 >eu = 0.02 (see in Govoni et 
al. discussion and comparison with other samples of radio galaxies). This is consistent with 
the idea that galaxies hosting BL Lac nuclei and/or radio sources have indistinguishable 
structures and have undergone a similar degree of interaction. 

In the microlensing scenario (Ostriker & Vietri 1985, 1990), what appears to be the 
host galaxy is actually an intervening galaxy whose stars are microlensing a background 
quasar. Since the alignment will not be perfect the galaxy will generally be off-center with 
respect to the point source, typically by ^ 0.5 arcsec (Merrifield 1992). To investigate this 
possibility we computed for each object the limit of off-centering, A = \(x, y)psF — (x, y)gai\, 
comparing the PSF location with that of the galaxy (derived by averaging of isophotal 
centers after subtraction of a scaled nuclear point source and exclusion of the central 
circle of radius 0.5 arcsec). With HST images the off-centering can be computed with 
unprecedented accuracy given the significantly better spatial resolution. The centers of 
nucleus and the host galaxy have a typical uncertainty of 0.2 and 0.4 pixels, respectively, so 
the net uncertainty is typically ~ 0.05 arcsec. 

None of the objects in our (low-redshift) sample show significant (> 2a) displacement of 
the point source with respect to the galaxy (values of A are given in Table 1); the observed 
displacements are always < 0.1 arcsec and therefore consistent with zero off-centering. 
Since the low-redshift BL Lacs were the prime candidates for microlensing according to 
Ostriker & Vietri (1985), our result strongly rules out the microlensing hypothesis for BL 
Lac objects in general. 

3.4. Isophotal Twisting 

For each host galaxy having ellipticity larger than 0.15 we have estimated the isophote 
twisting, A9, computing the difference of position angles over the whole range of observed 
surface brightness (see Table 1). For the others the small ellipticity introduces large errors 
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on the position angle so that it is not possible to get a reliable measurement of the isophotal 
twist. 

The largest observed twist is ~ 30° in 1853+671; this value is however uncertain 
because of the relatively small ellipticity and short exposure time. We found five other 
objects (corresponding to ~20% of the sample) with isophote twists larger than 15°. This is 
consistent (see Fig. 6) with what was found for a sample of 43 isolated (normal) ellipticals 
(Fasano & Bonoli 1989) and for 79 low-redshift radio galaxies (Govoni et al. 2000) using 
the same method to derive APA. 

A small amount of isophote twisting may be explained simply by the triaxiality of the 
galaxies, while larger twists may be due to the effects of tidal interaction with companion 
galaxies (Kormendy 1982). The lack of large twists in our sample is consistent with the 
overall smoothness (unperturbed shape) of the host galaxies. 

3.5. Isophotal Shapes: Disky and Boxy 

To check for small deviations from purely elliptical isophotes we analyzed the amplitude 
of the fourth cosine component, C 4 , of the Fourier fit to the isophotes (e.g., Bender & 
Saglia 1998). A significant positive value of C 4 corresponds to disky-shaped isophotes while 
a negative value indicates a boxy-shaped structure. 

We found that 80% of the sources show a C 4 amplitude (at the effective radius) 
smaller than 1% while no object exhibits |C 4 | larger than 3%. The distribution of C 4 , 
shown in Figure 7, is symmetric around zero (mean value < C 4 >bll= 0.03%, rms 
dispersion 0.7%) and since the size of the error is around 0.5% we conclude that in these 
galaxies there is no systematic deviation of isophotal shape from a simple ellipse. Similar 
distributions are observed in a sample of low-redshift radio galaxies using the same method 
(< C 4 >rg= 0.04% ±1.2%; Govoni et al. 2000) and also for radio galaxies in rich clusters 
(Ledlow and Owen 1995), as well as for normal ellipticals (Jorgensen et al. 1995). 

3.6. Structures in the Host Galaxies 

A first look at the raw images of BL Lacs as well as the surface photometry analysis 
above indicates that galaxies hosting BL Lacs are rather smooth and unperturbed 
(see images in Paper I). However, faint sub-structures could be hidden by the smooth 
contribution of the stellar component, so we subtracted a two-dimensional model of the 
galaxy using the best-fit parameters from the surface photometry analysis. This procedure 
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enhances faint structures such as companions, jet-like features, or any other high-contrast 
feature superposed onto the image of the galaxy (e.g., Faundez-Abans & Oliveira-Abans 
1998, and references therein). Clearly any region masked out during isophotal analysis 
(including the diffraction pattern of PSF spikes) will come out from this procedure. 
Moreover due to the sharpness of the nucleus and the insufficient sampling of the PC it 
was not possible to model the PSF adequately with this technique and some residuals 
(typically diffraction spikes) remain close to the nucleus. This has very little impact on our 
conclusions since structures very close to the nucleus are not considered. 

In Figure 8 we show three examples of the images of the BL Lacs before (left) and after 
(right) subtraction of the galaxy model. Panel (a) shows the optical jet of PKS 0521-36 
after subtraction of the galaxy-plus-point-source model. Panel (b) shows the image of BL 
Lac itself (2200+420) which, after subtraction of the model, does not exhibit significant 
features. Finally panel (c) shows an example (H 2356-309) of a companion at ~ 1 arcsec 
south-west of the nucleus. In all three subtracted images the emission from the diffraction 
spikes (not modeled) is clearly visible. 

The BL Lac objects in the present sample span the redshift range 0.03 £ z <^ 0.2. The 
usable field of view of the PC camera corresponds thus to projected linear size of ~ 25 
to ~ 135 kpc. At HST resolution therefore we are practically able to explore structures 
in the host galaxies as close as 0.2 - 0.9 kpc (depending on z) to the nucleus (we assume 
that within 5 pixels from the center the signal is strongly dominated by the nuclear point 
source) . 

In three cases (0521-36, 3C 371 and 2201+044) an optical jet is clearly visible (see 
Scarpa et al. 1999b for a detailed study of these objects). One object (0806+524) has 
arc-like emission at ~ 2 arcsec from the nucleus, which could be a shell around the galaxy, 
a feature not uncommon among ellipticals (see Scarpa et al. 1999a for details). No other 
faint structures (shells, bars or X features) are detected in the observed sources. 

HST observations of nearby radio galaxies have shown that dusty features are present 
in the central regions of most of the observed objects (Kleijn et al. 2000). These features 
often take the form of disks or lanes but sometimes the distribution is irregular. The 
observed dust pattern is generally confined within a region typically less than 1.5 kpc (in 
most cases the physical sizes of disks and lanes range from 0.3 to 1.5 kpc). Since BL Lacs 
are believed to be radio galaxies with the jet pointing toward us, we should observe similar 
dusty features also in our objects. However, two difficulties make this detection less likely 
in our sample objects. First, the bright nucleus tends to out-shine the central regions of the 
host galaxy. Second, the average distance of our objects is larger than that of comparison 
samples of radio galaxies observed with HST, making the angular size of dusty features 
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smaller. At z — 0.05, 1 kpc corresponds to ~ 0.75 arcsec; at this distance from the center 
the light from the nucleus is normally fainter than that from the host galaxy (see Fig. 1 in 
Paper I) so that these features could in principle be observable at the very lowest redshift 
but become progressively more difficult to detect at higher z. 

In order to explore this possibility we have therefore considered only the 9 objects 
at z < 0.1 and searched for dusty features in the central regions after subtraction of the 
bright point source from the original image. In only one object (1959+650), which is at 
redshift z = 0.048, is a significant dust lane observed at ~ 0.8 arcsec (corresponding to 
~ 1 kpc) north of the nucleus and roughly aligned with the major axis of the galaxy. This 
feature was also noted by Heidt et al. (1999) on ground-based images and is discussed 
in our previous work (Scarpa et al. 1999a) on peculiar objects observed in the snapshot 
survey. For the rest of the eight z < 0.1 sources no clear signature of dust is found. Note 
that four of these objects are at a redshift similar to or smaller than that of 1959+650. This 
lack of detection suggests that the relevance of dusty features may be different in BL Lacs 
and radio galaxies, with the caveat of the small statistics considered. Moreover, note that 
looking at features close to the nucleus requires not only high spatial resolution but also 
a PSF with sufficient sampling and faint wings. HST with the new Advanced Camera for 
Surveys High Resolution Camera offers two-times better sampling, and has a coronograph 
that depresses the PSF wings by factors of a few (and obscures the central 0.9 arcsec), so it 
could contribute to elucidating this point. 

4. Close Companions 

In a number of objects we noted the presence of faint companions around the target, 
a feature that seems also common to the HST images of quasars. We have systematically 
explored the near environments of the 30 low-redshift BL Lac sources with the aim 
of identifying potential companion objects. For three cases (0706+592, 1440+122 and 
2356-309) a very close (A < 1".2) compact companion is detected. These companions are 
located at projected distances of the order of 1-5 kpc from the nucleus. For five other 
BL Lac objects (see Table 2) diffuse or compact companions are found within a projected 
distance of £ 20 kpc (if they are at the same redshift as the BL Lac). 

In order to evaluate the statistical significance of the occurrence of these close 
companions we derived the expected number of objects within a given radius from the BL 
Lac based on the number counts observed in our PC frames well away from the BL Lac. 
The observed average density of objects with 19.5 < < 23.5 mag is 9 ± 5 arcmin" 2 , in 
good agreement with the observed galaxy density in several high galactic latitude fields 
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(Metcalfe et al. 1991). 

At these faint magnitudes the counts at high galactic latitudes are dominated by 
galaxies. The average galactic latitude of our fields is |6| = 45 deg but two fields (2200+420, 
2344+514) are at low galactic latitude (|6| ~ 10 dcg ) and one (1514-241) is at \b\ ~ 27 dcg but 
close to the direction toward the galactic center. The observed average density of stellar 
(unresolved) objects including all PC frames is 3.2 ± 4.2 arcmin" 2 , while excluding the 
three fields mentioned above yields a star density of 2.1 ± 2.6 arcmin~ 2 . This is consistent 
with other estimates of the stellar density at high galactic latitude (Hintzen, Romanishin & 
Valdes 1991; Metcalfe et al. 1991). 

Assuming the observed average density of faint (19.5 < m R < 23.5 mag) objects, 
9 arcmin -2 , around 30 targets we would expect to find 0.06 companions at r <0'.'5, while 
one (a compact companion at 0.3 arcsec from 1440+122) is observed (Poisson probability 
P=0.06). Within 1'.'5 we would expect to observe 0.6 objects while 3 are detected (P=0.02). 
At larger radii (r <5") we observe 11 companions when 6 are expected (P=0.02). To test 
the reality of this apparent excess, we counted the number of detected objects within these 
same radii for random positions in each PC frame. In those cases we found 0, 1, and 7 
companions, respectively, in agreement with the expected numbers of faint field objects. 

Although the numbers are small and consequently the statistics are not very good, we 
conclude from our analysis that there is some indication that an association between BL 
Lacs and companions is not dominated by chance alignments. 

The identity or role of these companions is not yet clear, however. In other BL Lac 
fields, optical spectroscopy of the companion has sometimes proved the physical association 
with the active galaxy (Falomo 1996); in other cases, they have proved to be stars or 
foreground or background galaxies, fortuitously aligned with the BL Lac (e.g., Stickel et al. 
1991). For most cases here, spectra of the companion objects are not available and only 
statistical arguments suggest some are likely associated with the BL Lac object. 

Some of the close companions are resolved faint galaxies with absolute magnitudes 
(assuming they are at the distance of the BL Lac object) in the range —21 < m R < — 19 mag 
(Table 2). This is consistent with previous evidence that BL Lac objects inhabit 
environments with higher-than-average galaxy densities, as in groups or poor clusters 
(Pesce, Falomo & Treves 1995; Wurtz et al. 1996). In the two cases with more than one 
companion, 1229+645 and 1440+122, the two companion objects are clearly different from 
one another, one being resolved and the other unresolved, 

We note that a significant number of the observed companions are compact and 
appear point-like even at HST resolution. Based on the counts of Galactic stars, we would 
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expect only 1 or 2 chance alignments with stars while 6 compact companions are observed 
(Table 2). The excess of compact companions is in fact more significant than the excess of 
resolved companions (galaxies), especially since none of the compact sources belong to the 
low galactic latitude fields where the star density is higher than the average value. 

The nature of these unresolved companions is unclear. If they were at the redshift of 
the BL Lac object, their magnitudes would be in the range —18 < m R < —15 mag, too 
luminous for globular clusters. These magnitudes are consistent with dwarf ellipticals but 
the observed sources are probably too compact, unless they contain weak active nuclei (but 
they would have to be very weak indeed). One can ask whether some of these compact 
sources could be supernovae of type la. The SN la rate in early type galaxies is estimated 
to be about 0.1-0.2 SN/100yr/10 10 Lg (Cappellaro et al. 1999). This means that in a 
massive galaxy of M R ~ -23.7 mag we expect to have about 1-2 SN/lOOyr. Because the 
SN la at maximum reach M R ss -19.4 mag and our limiting magnitude for stellar objects 
is rriR ~ 24.5 mag (it is m R ~ 24 mag for extended objects), at z ~ 0.1 we should be able 
to detect such SN for about 100 days. Therefore we should find 0.2 SN in our 30 galaxies 
observed. Unless the SN rate in these active galaxies is significantly higher than the average 
for normal ellipticals, it is highly unlikely that the compact companions are SN. 

We believe the nature of all these companions, resolved or unresolved, will remain 
uncertain until spectroscopy is performed. Whatever their identity, there is little evidence 
of tidal interaction in the BL Lac host galaxy (appreciable isophote twisting, non-concentric 
isophotes, etc.), which suggests the companions are not strongly perturbing it. Gravitational 
interactions might have had a greater role in the past in the formation and fueling of the 
active nucleus, and what we see now could be a leftover of an earlier (close) interaction. 
One interesting possibility suggested by numerical simulations (Bekki 1999) is that the 
close companions are the product of a past major merger between gas-rich galaxies and 
also provide fuel for the nuclear activity but are not necessarily linked to the formation of a 
massive black hole in the nucleus, or to the sustenance of an active relativistic jet. 

5. Summary and Conclusions 

We have presented detailed morphological analysis of HST images for 30 galaxies 
hosting BL Lac sources at redshift z ^ 0.2. The galaxies were investigated with an order of 
magnitude better spatial resolution than it is possible from the ground. 

Our analysis indicates that in spite of the presence of the active nucleus the host 
galaxy appears in most cases to be a completely "normal" elliptical. Both the ellipticity 
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and the isophotal shape distributions are similar to those for radio galaxies and radio-quiet 
ellipticals. This suggests that tidal interactions are very infrequent or are short-lived with 
respect to the nuclear activity time scale. 

We find no indication of displacement and/or off-centering of the galaxy isophotes with 
respect to the nucleus, meaning the unresolved nuclear source truly sits in the center of the 
galaxy. This rules out the microlensing hypothesis for BL Lacs, which predicts frequent 
off-centering of the nucleus. This does not exclude the possibility that some objects could 
be lensed (see Scarpa et al. 1999a), but it cannot be a widespread explanation of the BL 
Lac phenomenon. 

An interesting result is the finding of excess close companions in several cases. A high 
incidence of companions of BL Lacs was first suggested by (Falomo, Melnick & Tanzi 1991) 
from sub-arcsec resolution ground-based imaging, and was later confirmed by Pesce, Falomo 
and Treves (1995), Falomo (1996) and Heidt et al. (1999). Close companions seem also to 
be common around quasars imaged by HST (Bahcall et al. 1997; Disney et al. 1995). 
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Structural Properties of Low-Redshift Host Galaxies 
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(°) EUipticity of the host galaxy at the effective radius. 

Typical error ± 0.02 far from the nucleus. 
^ Fourier coefficient 100*6*4 describing the isophote 

shape (see text); typical uncertainty is 0.5. 
Relative displacement of isophotes (dimensionless, 

see definition in § 3.3); typical uncertainty is ± 0.01. 
W Off-centering (arcsec) of the nucleus with respect to the host 
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Table 2 

Properties of Close Environments of BL Lac Objects 
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(a) Absolute magnitude assuming the same redshift as the BL Lac object. 

(b) Projected distance from the nucleus (in arcsec and kpc) assuming 
assuming the same redshift as the BL Lac object. 

(c) Position angle of the feature. 
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Fig. 1. — Distribution of host galaxy absolute magnitudes for the full HST snapshot sample 
(Paper II). Hatched area: The low-redshift subsample discussed in this paper. 
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Fig. 2. — Isophotes derived from the surface photometry analysis (dotted lines) and fitted 
ellipses (solid lines) for the BL Lac object 1426+428. 
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Fig. 3. — The radial surface brightness distribution of the BL Lac 1426+428 derived from 
the two-dimensional analysis (filled circles) agrees very well with the azimuthally averaged 
one-dimensional profile (solid line). The inset shows the brightness profile (filled circles) 
compared with the de Vaucouleurs fit (solid line). The excess at small radii is due to the 
presence of the bright nucleus (see Paper I for details). 
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Fig. 4. — Comparison of the distribution of ellipticities for the host galaxies of BL Lac objects 
(top) with those of normal ellipticals (middle) from a sample of nearby galaxies (Fasano & 
Vio 1991) and low-redshift radio galaxies (bottom; Govoni et al. 2000). 




Fig. 5. — Comparison of the isophote displacement, 5 (see text), for host galaxies of BL Lacs 
(top; present work) and low-redshift radio galaxies (bottom; Govoni et al. 2000). 
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Fig. 6. — Comparison of the maximum isophotal twisting, A6 max , for BL Lac host galaxies 
(top), 43 non-radio elliptical galaxies (middle; Fasano & Bonoli 1989), and low-redshift radio 
galaxies (bottom; Govoni et al. 2000). 
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Fig. 7. — Comparison of the C4 parameter (see text) for BL Lac host galaxies {top), normal 
ellipticals and SO galaxies (middle; Jorgensen et al. 1995), and low-redshift radio galaxies 
(bottom; Govoni et al. 2000). The distributions are similar, and consistent with little or no 
diskiness or boxiness. 
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Fig. 8. — Example images of low-redshift HST-observed BL Lac objects (left panels) 
compared to galaxy-model-subtracted images (right panels). Both panels are displayed using 
the same grey-scale. Top: PKS 0521-365 and its optical jet and bright knot ~ 2 arcsec 
from the nucleus. Middle: BL Lac (2200+420). No significant residuals are detected after 
subtraction of the model (point source plus elliptical galaxy, convolved with the PSF). 
Bottom: H 2356-309. A compact companion is detected 1.2 arcsec from the nucleus of 
the galaxy. The cross-shaped residuals are due to the diffraction pattern of the HST PSF, 
which is not fully described by the two-dimensional modeling. Field orientation (arrow mark 
North) and scale are given for each image. 
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